Numerical Methods of Applied Mathematics I

IEEE-standard, iterative and direct linear system solution methods, eigendecomposition and model-order reduction, fast Fourier transforms, multigrid, wavelets and other multiresolution methods, matrix sparsification. Nonlinear root finding (Newton's method). Numerical interpolation and extrapolation. Quadrature

A scalable FETI-DP algorithm for a semi-coercive variational inequality

We develop an optimal algorithm for the numerical solution of semi-coercive variational inequalities by combining dual-primal FETI algorithms with recent results for bound and equality constrained quadratic programming problems. The discretized version of the model problem, obtained by using the FETI-DP methodology, is reduced by the duality theory of convex optimization to a quadratic programming problem with bound and equality constraints, which is solved by a new algorithm with a known rate of convergence given in terms of the spectral condition number of the quadratic problem. We present convergence bounds that guarantee the scalability of the algorithm. These results are confirmed by numerical experiments

A scalable FETI-DP algorithm for a coercive variational inequality

We develop an optimal algorithm for the numerical solution of coercive variational inequalities, by combining FETI algorithms of dual-primal type with recent results for bound constrained quadratic programming problems. The discretized version of the model problem, obtained by using the FETI-DP methodology, is reduced by the duality theory of convex optimization to a quadratic programming problem with bound constraints. The resulting problem is solved by a new algorithm with a known rate of convergence given in terms of the spectral condition number of the quadratic problem. We present convergence bounds that guarantee the scalability of the algorithm. These results are confirmed by numerical experiments